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1. Introduction 


Most Computational Fluid Dynamic (CFD) codes used in practical applications are based on the 
Reynolds-Averaged Navier-Stokes (RANS) set of equations. These equations employ many ad hoc 
assumptions and empirical coefficients in the modeling of the turbulent viscous-fluxes and the chemical 
specie-production source-terms. Therefore, these models have to be validated with experimental data 
before being applied with any confidence to a given class of flows. The types of flow of interest for the 
present work arc those found in scramjet applications. 

A supersonic-combustion scramjet model known as SCFIOLAR has been tested at NASA Langley’s 
Direct-Connect Supersonic Combustion Test Facility (DCSCTF) ’ . This experiment was designed to 
provide a test case for validation of CFD codes, involving supersonic injection, mixing and combustion 
in a duct. The model has a simple geometry, and large regions of subsonic/recirculating are avoided; the 
enthalpy of the test gas corresponds to a Mach-7 flight condition (approximately 1200 K at combustor 
entrance). During the design of the experiment it was believed that the resulting flow would be mixing- 
limited (i.e., with chemical reaction towards equilibrium being much faster than the mixing process); as 
will be seen later, this turned out not to be the case. Some numerical simulations of this experiment have 
already been reported by other authors . 

The present work documents some preliminary results in the numerical simulation of the 
SCFIOLAR scramjet using VULCAN, a RANS computer code developed at NASA Langley 4 . The pur- 
pose of the simulation is to assess the effects of the modeling assumptions and parameters in the flow- 
field predictions. A brief description of the experiment and summary of the data is presented first. It is 
followed by the numerical analysis of the facility nozzle, which will provide the necessary inflow condi- 
tions for the simulation of the direct-connect combustor. The analysis of the scramjet proper begins with 
the definition of a baseline configuration for the numerical modeling. This configuration provides a ref- 
erence solution which is used for a qualitative analysis of the flowfield, grid convergence study, and 
comparison with the data. The main objective of the investigation is to quantify the effects on this solu- 
tion due to changes in some of the numerical parameters. This report concentrates on turbulence param- 
eters such as numerical models, Schmidt and Prandtl numbers, and inflow turbulence. Future work will 
include the effects of chemistry models, inflow conditions and the presence of radicals. 



2. Experiment 


The SCHOLAR experiment has been extensively described in references 1 and 2. Only a brief 
description of the experimental set-up and a summary of the data will be presented here. 

2.1 Description of the Facility 

In the DCSCTF facility, vitiated air is produced at high enthalpy in the heater by combustion of 
hydrogen with premixed oxygen and air; the nominal enthalpy corresponds to Mach 7 flight. This air is 
accelerated through a water-cooled Mach-2 nozzle before entering the test section (operating conditions 
will be given later in the report). The complete layout and dimensions of the combustor model are 
shown in figure 1 (a). The duct consist of a constant-area isolator, which ends in a small step on the top 
wall. Downstream of the step there is another short constant-area section. This section is followed by a 
3° divergence on the top wall, which extends all the way to the exit of the model. The span of the model 
is constant throughout its entire length. The injection nozzle provides hydrogen fuel at a 30° angle [fig- 
ure 1 (b)]. This nozzle was designed to produce a Mach 2.5, one-dimensional flow at its exit (injection 
conditions and equivalence ratio to be given later). Also shown in the figure are 5 pilot injectors that 
were used in some runs to pilot the main injector; these injectors were not modeled in the present calcu- 
lations. 

The model was built in two sections: the upstream section (containing the isolator and injector) was 
made of copper, while the downstream section (containing most of the divergent section) was made of 
steel. The conductivity of these materials and the thickness of the walls allowed the duct to be run with- 
out cooling. There are 7 transverse slots in the model to allow for measurement of static temperatures 
fields by CARS beams; the locations shown in figure 1 (numbered 1, 3, 5, 6, 7) are the ones used in this 
study. The model itself is instrumented with both pressure taps and wall-temperature probes. Pressure 
taps are located at the centerline of the bottom wall, on the top wall (close to the sidewall for the copper 
section and centerline for the steel section), and on the sidewall (midpoints of the steel section). Ther- 
mocouples are located on the top wall, six for each section of the duct. 

2.2 Summary of Data 

Wall-pressure distributions for the pressure taps at the bottom centerline are shown in Figure 2. 
Data are shown at 10 s and 24 s; note that the heater is started at 1 s and fuel injection commences at 6.4 
s. Comparison between measurements at these two times reveals only small differences. Also shown is 
the average over the total time run; this average will be used throughout the present work for compari- 
son with numerical results. Pressures vary widely in the upstream region, probably due to the shock- 
wave system created by the injector and step; a small discontinuity in cross-section between nozzle and 
combustor, due to a stainless-steel flange, may also account for the pressure variations. Pressure gener- 
ally falls moving downstream due to the divergence in the duct, until about 0.5 m. At that location it 
rises rapidly to a peak at about 0.75 m, and then falls slowly as the gas flows towards the exit. Very little 
combustion appears to take place upstream of 0.5 m, but by 0.75 m most of the mixed fuel seems to 
have burned; further mixing and combustion is expected to occur downstream. 

Since the duct is uncooled, surface temperatures vary greatly during the course of the run; represen- 
tative temperature histories and fits (top wall, unpiloted) are presented in figure 3. In the copper section, 
temperature is typically about 350 K at the beginning of the fuel injection and rises to as high as 600 K 
at the end. For the steel section, those numbers are typically 450 K and 950 K, respectively. 
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Figure 4 contains three-dimensional cutaway views of the duct showing contour-plots of the fitted 
CARS temperature data; the data are all fit to cosine-series bivariate functions. Flow-direction is from 
top-left to bottom-right; main fuel injector is on the top wall between planes 1 and 3. Planes 3 and 5 
show a region of low temperature (between 250 and 550 K) at the center, which is the injected fuel- 
plume; there is no evidence of combustion between these two planes. At planes 6 and 7 temperatures 
have risen abruptly, suggesting near complete combustion. The hot regions at the top and bottom walls 
are probably near stoichiometric, with the cooler region at the center likely to be fuel-rich. Injected fuel 
has probably not penetrated to the sidewalls. There is some asymmetry in the data, probably due to 
inflow perturbations and the data fit. 
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3. Analysis of the Facility Nozzle 


As mentioned, a two-dimensional, converging-diverging, water-cooled nozzle was used to acceler- 
ate the vitiated flow to the nominal test conditions at the combustor inlet. The conditions at the nozzle 
exit were probed with a pitot rake to map the pitot pressure (po 2 )- These measurements were reported in 
reference 5, together with comparisons to CFD calculations performed with the VULCAN code. The 
facility nozzle was modeled in the present work to provide inflow conditions to the duct calculations. 

3.1 Solution Procedure 

The three-dimensional grid used for the nozzle simulation is shown in figure 5; its dimensions are 
129 x 129 x 145 (axial x vertical x lateral). The half-width and full height of the nozzle was modeled, 
the latter because of slight asymmetries between top and bottom (curved) walls. The grid was subdi- 
vided into 32 blocks, to take advantage of the Message Passing Interface (MPI) capabilities of VUL- 
CAN. It had a wall-spacing of 5 x 10" 5 m, resulting in a maximum y + of 65. 

The inlet boundary was modeled as subsonic inflow, with a total pressure and temperature of 
767250 Pa and 1828 K, respectively. The inflow turbulence intensity x and turbulent-to-laminar viscos- 
ity ratio P [/p[ were set at 25% and 600, respectively (which result in a turbulence length-scale of about 
1 mm at inflow conditions). All solid (water-cooled) walls were modeled as no-slip isothermal bound- 
aries; following reference 5, the wall-temperature was assumed to be 500 K. The exit boundary was set 
as supersonic extrapolation. Finally, the vertical centerplane was modeled as a symmetry boundary. As 
mentioned, the inflow to the nozzle was vitiated air. One-dimensional calculations assuming equilib- 
rium and frozen compositions differed little in major species, temperature and pressure 1 . Because of this 
and the range of temperatures under consideration, the gas was modeled as calorically perfect in order to 
reduce computational expense. Following the one-dimensional results, the molecular weight and spe- 
cific-heat ratio y were assumed to be 25.9098 and 1.28909, respectively. 

The inviscid fluxes were calculated using Edwards’ low-dissipation flux-split scheme 6 , with third- 
order MUSCL interpolation. Time-integration was done using the implicit diagonalized approximate- 

n 

factorization scheme , the entire domain was run fully elliptic. Wilcox’s 1998 k-co model and wall- 

o 

matching functions were used for modeling turbulent viscous stresses and heat-transfer; the turbulent 
Prandtl number (Pr T ) was assumed to be 0.90. Calculations were done on an SGI Origin 2000, using 32 
processors. Grid sequencing with three levels was used for convergence acceleration, and also as a mea- 
sure of grid convergence, as will be seen in the results. Convergence was considered to have been 
achieved when the residual dropped four orders of magnitude, which occurred after 1500 iterations in 
the fine sequence. 

3.2. Results 

Figure 6 shows the nozzle centerplane Mach contours. The effects of the wall asymmetry are just 
apparent in the contours. The calculated nozzle-exit pitot-profiles are presented in figure 7 for the verti- 
cal centerplane, and in figure 8 for the horizontal centerplane. Both show adequate agreement with the 
experimental data, although the latter one shows more asymmetry and thicker boundary-layers (espe- 
cially in the horizontal centerplane). It should be noted the data was obtained at slightly different operat- 
ing conditions than those of the SCFIOLAR test used here. Also, the pitot rake was located 
approximately 2 mm downstream of the nozzle exit (with the duct removed); the numerical results were 
taken at the exit itself. Finally, measurements may be affected by the probe diameter not being small 
enough with respect to the boundary-layer thickness. 
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A detail of figure 7 is shown in figure 9 for the medium- and fine-sequence calculations, together 
with 1% error bars on the data. The fine-sequence solution falls within the error bars for most of the 
points, although the data shows a greater dip near the center. There are some differences between 
medium and fine calculations, but for the most part these differences are of the order of the experimental 
error; therefore, the solution can be considered to be grid-converged. 

The values of the inflow turbulence-parameters quoted before were chosen in order to provide 
scramjet-inlet turbulence-intensity and viscosity -ratio of the order of 1% and 100, respectively (giving a 
turbulence length-scale of the order of 1 mm at exit conditions). There is no way to quantify the turbu- 
lence velocity and length scales of the experiment. To determine the effects of these parameters on the 
scramjet simulation, several combinations were chosen. It was found difficult to predict a priori which 
inflow values would give desired values at the exit, but after some trial-and-error the inflow parameters 
shown in table 1 were found to give adequate values at the exit potential flow (also given in the table). 
In particular, the values of the inflow turbulence parameters given earlier provide exit conditions with 
2.5% turbulence intensity and 300 viscosity-ratio. The effects on the nozzle-exit vertical-profiles of the 
different combinations of turbulent parameters given in the table can be seen in figure 10. In general, 
higher values of x and p^/p^ appear to give better agreement with the nozzle experimental data. The 
effects of these parameters on the flow within the combustor proper will be shown in the next section. 
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4. Scramjet Baseline Solution 


4.1 Solution Procedure 

A reference solution was sought that would provide a first approximation to the problem, and from 
which parametric studies could be done. This solution, and the procedure and conditions used to obtain 
it, will hereafter be referred to as “baseline”. Unless otherwise noted, this numerical configuration have 
been used for all calculations reported here. 

The half-width combustor was discretized with a grid that had approximately 3.6 million control 
volumes (CVs), uniformly distributed among 80 blocks (again, to take advantage VULCAN’s MPI 
capability). A detail of the grid in the neighborhood of the injector is shown in figure 1 1. To reduce grid 
requirements, non-C(O) block-boundaries were extensively used (as may be seen from the grid disconti- 
nuities in the figure). The vertical grid-dimensions varied between 61 and 85, and the lateral between 61 
and 81. The grid spacing at the wall varied between 1.5 x I (T 4 m and 3.0 x 10“ 4 m. For the most part, the 
resulting values of y + are below 100, which is just about right for the use of wall functions. 

For inflow conditions at the combustor inlet, the profiles obtained from the nozzle calculation were 
imposed on the boundary. To allow for multi-species calculations, species mass-fractions were added. 
The vitiated composition from the one-dimensional calculations previously mentioned included about 
two dozen species. For the chemistry models used in the present calculations (which could not account 
for most of the vitiation species), the inflow was assumed to contain only N 2 , 0 2 and F1 2 0 with uniform 
mass-fractions 0.5638, 0.2321 and 0.2041, respectively 3 ; these were obtained from the one-dimensional 
analysis using the same mole fractions for 0 2 and F1 2 0 (the mole fraction of N 2 was slightly adjusted to 
give a total mole-fraction of 1.0). The inflow conditions at the inlet of the hydrogen injector were 
assumed to be subsonic with total pressure and temperature of 3.44 MPa and 302 K, respectively. These 
conditions resulted in an equivalence ratio of approximately 1 .0. The solid walls were modeled as no- 
slip isothermal, with wall temperatures set at 400 K and 500 K for the copper and steel sections, respec- 
tively (based on the wall-temperatures shown in figure 3). The exit boundary was modeled as super- 
sonic extrapolation, and the vertical centerplane was modeled as a symmetry boundary. 

The inviscid space and time integration schemes were the same as for the nozzle calculations, 
unless otherwise noted. The burner was ran fully-elliptic, to account for possible subsonic/recirculation 
regions; some calculations were done using a semi-elliptic approach to be described below. The flow 
was modeled as a mixture of thermally-perfect gases. The chemistry model used was Drummond’s 9- 
specie/ 18 -reaction (9x18) mechanism 9 ; the ignition-region feature in VULCAN was needed to ignite 
this model at the present conditions. Thermodynamic properties of the chemical species were modeled 
with a one-curve polynomial fit included in VULCAN, and based on the data from reference 10. The 
Wilcox k-co turbulence model and wall functions were again used. They were supplemented with the 
Wilcox compressibility correction, used to reduce spreading rates (or turbulence) as a function of the 

o 

local turbulence Mach number . Not included was the Pope correction for the round-jet/plane-jet anom- 
aly, which is standard in Wilcox’s model. The turbulence Prandtl (Pr T ) and Schmidt (Sc T ) numbers 
were set at 1 .0. 

A three-level grid sequencing was again used for convergence acceleration and to estimate grid con- 
vergence. Computations were done on 80 500-MFIz processors of an SGI Origin 2000 computer; the 


a The effects of radicals such as OH present in the inflow will be the subject of future work. 
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ideal parallel speed-up obtained with the present grid was 78.3. Also used were 16 DEC Alpha 667 
MHz nodes in a BEOWULF cluster. 

4.2. Results 

A total of 10,000 iterations were done on the medium sequence, and 45,000 on the fine sequence. 
For both sequences, the residual would drop about two orders of magnitude before undergoing slow 
oscillations. For the medium-grid solution, convergence was assumed when no variation would be 
observed in either the overall flowfield or the wall-pressure distributions; the solution was mostly con- 
verged after 2,000-3,000 iterations. For the fine-grid solution, a minimum of 35,000 were needed before 
the solution would show no further changes. 

Figure 12 shows the centerplane Mach contours for the baseline solution (note that in this and sub- 
sequent centerplane plots the vertical scale is twice the horizontal for better resolution). The jet pene- 
trates at approximately 30° for a certain distance, but it is immediately turned in the direction of the 
main flow while continuing to penetrate further towards the opposite wall. The expansion after the step, 
and compression from the jet are clearly seen. A region of slower flow is apparent downstream of x ~ 
0.50 m. 

The water mass-fraction contours (figure 13) suggests that ignition takes place between 0.45 and 0.5 
m. This is also shown by the temperature contour plots (figure 14). The CARS measurements samples 
only indicate that heat-release occurs after 0.4 m. The fuel core (indicated by the low-temperature 
regions) seems to penetrate further in the numerical simulation. The mixing and combustion efficiency 
distributions are presented in figure 15; they are defined 11 as mixed- fuel over total- fuel, and reacted 
(water)-fuel over total-fuel, respectively, and corrected for vitiation . Initially, the combustion lags 
with respect to the mixing, but it raises rapidly after 0.5 m, and by 0.8 m most of the mixed fuel appeal's 
to have reacted. Overall, just over 60% of the fuel seems to mix and react with the inflow air. 

These observations are confirmed by the analysis of the bottom-wall centerline pressure distribution 
(figure 16). The experimental data sampling can only locate the ignition delay somewhere between 0.45 
and 0.57 m. CFD predicts the flame to begin near the upstream limit; this is likely to be an underpredic- 
tion of the delay. The pressure levels from the medium- and fine-sequence solutions are compared in 
figure 17. At the inlet, the medium- solution levels are slightly higher than the fine-solution; this has 
been observed in other VULCAN calculations , and may be due to post-processing of the medium- 
sequence solution. The main difference between both solutions is in the ignition delay, which is consid- 
erably greater for the medium solution (and closer to the data). This is counter-intuitive, since the 
greater numerical diffusion of the coarser grid should result in greater mixing, and therefore less igni- 
tion delay (see below). In any case, these results show that grid convergence at the very least has not 
been proven on the current grid. 

This preliminary solution indicates that the baseline model likely undeipredicts the ignition delay, 
or time for a fuel-air mixture to form a flame in the absence of ignition sources 14 . This delay is the com- 
bination of a physical delay (depending on fuel mixing) and a chemical delay (time for flame production 
after the formation of a homogeneous fuel-air mixture). In the numerical solution, the former mostly 
depends on the turbulence model and the latter on the chemistry model. Therefore, the underprediction 
of the ignition delay may be due to an overprediction of the mixing by the turbulence model, an under- 
prediction of the chemical delay by the chemistry model, or a combination of both. Another deficiency 
is the overprediction of the penetration of the fuel plume; it is likely due to the turbulence model, 
although this is subject to further verification. 
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Some of the following calculations were obtained with a “semi-elliptic” integration approach, in 
order to reduce computational requirements. In this approach, the first 48 upstream blocks were solved 
elliptically, with a first-order extrapolation boundary-condition applied at the downstream end of the 
resulting region (at x ~ 0.35 m). Convergence on this region was fairly rapid, owing to the absence of 
significant reaction. The remaining 32 blocks were then solved, also elliptically, for about 15,000- 
20,000 iterations. The results for the baseline condition are shown in figures 18 and 19. The solution is 
almost indistinguishable from the fully-elliptic calculation. Since most of the computation is spent in 
the second (reacting) region, this approach reduces the computational effort by almost 60%. The calcu- 
lations performed with this approach will be explicitly noted. 

As explained in the Introduction, for the present work the effects in the turbulence parameters were 
studied. Two additional turbulence models were tested: Menter’s k-co model 15 , and an explicit algebraic 
stress model extension (EASM) of the Wilcox model based on Gatski and Speziale’s model ; these 
calculations were done using the semi-elliptic integration described above. Menter’s model combines 
the Wilcox k-co model near solid walls, with the standard k-e model everywhere else (particularly in 
shear-dominated flow-regions). The results for this model (figures 20 and 21) show that it slightly over- 
predicts ignition delay, and has less fuel-plume penetration than the Wilcox model. On the other hand, it 
considerably overpredicts the pressure-rise in the reaction region. The EASM model is essentially an 
extension of the Boussinesq approximation, providing a non-linear constitutive relation between turbu- 
lent stresses and strain-rates. They account for the anisotropy of the normal Reynolds stresses, which in 
turn allows for the simulation of stress-induced motions. These motions are of importance in corner- 
flow regions, were they may delay the onset of separation due to adverse pressure gradients. For the 
present application, it can be seen (figures 22 and 23) that the EASM model provides an acceptable 
approximation to the ignition delay and pressure-rise, but at the expense of overpredicting the fuel- 
plume penetration. 

The effects of lowering Sc T and Pr T by half from their baseline values are shown in figures 24 to 27. 
Lowering Sc T (figures 24 and 25) appears to have little impact on fuel-plume penetration, but the igni- 
tion-delay and pressure -rise predictions are unacceptable. On the other hand, a lower Pr T has little effect 
on pressure -rise and ignition-delay while having a somewhat lower penetration of the fuel-plume. The 
values Pr T = Sc T = 0.50 are usually recommended for shear-dominated flows . No attempt was made 
to test values higher than the baseline ones, since it was considered that these were already high enough. 

The effects of the different combinations of inflow turbulence- parameters of table 1 are given in fig- 
ures 28 to 31. The high-turbulence inflow (figures 28 and 29), which appeared to give better agreement 
with the nozzle data, is shown here to give a somewhat better approximation to jet-penetration, at the 
expense of the agreement with ignition delay and pressure-rise. The low-turbulence inflow (figures 30 
and 31) seems to overpredict ignition delay and penetration, and somewhat underpredict pressure -rise. 
In brief, a change of approximately 50% in x and a factor of 3 in p T /p L with respect to baseline values 
can give drastically different solutions. These parameters are not usually estimated in experiments, and 
in any case can only be done so in an order-of-magnitude basis. 

One common thread through all these calculations is the apparent inability of current turbulence- 
modeling to accurately predict ignition delay, combustion pressure-rise and fuel-plume penetration 
simultaneously. So far, only two of the three parameters at best could be adequately captured. 


b It should be noted that this model is still under development within VULCAN. 
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5. Conclusions 


The SCHOLAR scramjet is the subject of an ongoing numerical investigation, with the puipose of 
validating the VULCAN code and its turbulence and chemistry models. The facility nozzle was mod- 
eled first, in order to provide adequate inflow conditions to the combustor. Results for the latter show 
that CFD could predict its qualitative behavior, particularly reaction-limited heat-release due to ignition 
delay. From a quantitative point of view, at least two problems were identified: the inability to prove 
grid-convergence for the grid currently used, and the difficulty of accurately predicting all three param- 
eters relevant for the present application: ignition delay, combustion pressure-rise, and fuel-plume pen- 
etration. 

Regarding grid convergence, this could not be proved using simple grid-sequencing. Additionally, 
the trends between the sequences were not what would be expected in terms of numerical diffusion. One 
possible cause for this behavior could be in the extensive use of non-C(O) zonal boundaries, both normal 
and parallel to the flow. The latter ones could adversely affect the numerical flowfield by introducing 
extended areas of numerical diffusion. 

The second problem is associated with the accuracy of turbulence and chemistry models. Only the 
first ones were investigated here, and the numerical solution was shown to be sensitive to them. Differ- 
ent turbulence models gave qualitatively dissimilar solutions. Furthermore, a relatively high value of 1.0 
had to be used for the turbulent Schmidt number, if the ignition delay was to be captured at all. Another 
issue is that of inflow turbulence-parameters. These are not usually determined in experiments, and 
appear to have a considerable impact in numerical solutions. It is not immediately obvious what combi- 
nation, if any, of modeling parameters would allow to simultaneously predict the three flow-parameters 
mentioned above with any acceptable accuracy. 

Future work will initially concentrate on obtaining a grid-converged solution. A new grid is under 
development in which the discontinuous boundaries have been restricted to a limited number normal to 
the flow; the computational price to pay is an almost doubling of the total number of CVs. Based on the 
results from this grid, it will be determined whether the numerical results obtained so far are valid or 
need to be recalculated. The next step will be to include the effects of chemistry parameters as described 
in the Introduction. The overall objective will be to determine the sensibility of the solution to the dif- 
ferent parameters present in most numerical models, and perhaps establish an “envelope” in some solu- 
tion space were the “correct” answer may be expected. 
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Tables 


Table 1 . Nozzle Inflow and Outflow Turbulence Parameters 


Description* 

Inlet 

Exit f 

X 

Bt/Bl 

X 


Baseline 

0.25 

600 

0.025 

300 

High turbulence 

0.10 

600 

0.035 

850 

Low turbulence 

0.25 

300 

0.015 

100 


Based on exit conditions 

'Values are approximate potential-core averages. 
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Figures 


(a) 




Figure 1. Layout of the combustor model: 

(a) Nozzle and complete combustor, (b) Detail in the region of the injector. 
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Figure 3. Top-wall temperature history. 
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Figure 4. Cutaway planes with contours of mean temperature 
(dimensions in meters). 
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Figure 5. Three-dimensional nozzle grid (every other grid-line shown). 
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Figure 6. Nozzle centerplane Mach contours. 
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Figure 9. Nozzle vertical-centerplane pitot-pressure profiles: grid sequencing. 



Figure 10. Nozzle vertical-centerplane pitot-pressure profiles: effects of inflow turbulence. 
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Figure 12. Centerplane Mach contours - baseline solution. 
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Figure 13. Centerplane water contours - baseline solution. 
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Numerical 

Figure 14. Temperature contours - baseline solution. 



Figure 15. Efficiencies distributions - baseline solution. 
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Figure 16. Bottom-wall centerline pressures- baseline solution. 



Figure 17. Effect of grid sequencing on wall-pressures - baseline solution. 
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Figure 18. Water contours - Semi-elliptic integration. 
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Figure 21. Wall-pressure distribution - Menter’s turbulence model. 
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Figure 23. Wall-pressure distribution - EASM turbulence model. 
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Figure 24. Temperature contours - Sc = 0.50. 



Figure 25. Wall-pressure distribution - Sc = 0.50. 
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Figure 26. Temperature contours - Pr = 0.50. 



Figure 27. Wall-pressure distribution - Pr = 0.50. 
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Figure 28. Temperature contours - High inflow-turbulence. 



Figure 29. Wall-pressure distribution - High inflow-turbulence. 
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Figure 30. Temperature contours - Low inflow-turbulence. 
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